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LONG-TERM  GOAL 

To  develop  an  inversion  procedure  which  will  estimate  the  bottom  acoustic  properties  with  sufficient 
resolution  to  enable  prediction  of  acoustic  field  over  a  range  of  frequencies. 

OBJECTIVES 

Investigate  the  usefulness  of  wavelets  in  determining  the  bottom  compressional  wave  speed  profile  at 
different  scales. 

BACKGROUND 


The  wavelet  y/x  u  (l  )  represents  a  family  of  functions 
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where  X  is  the  parameter  that 


controls  the  dilation  or  contraction  of  the  function  and  u  the  translation  of  the  function.  In  practical 
applications  the  scale  parameter  X  and  translation  parameter  u  are  discretized.  A  discrete  representation 


of  the  wavelet  is  y/m  n  (l)  = 
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where  m  is  the  scale  parameter  and  n  is  the  position 


parameter.  For  certain  class  of  wavelets  the  function  i//m  n  (/  )  are  orthonormal  to  their  dilates  and 


translates.  These  wavelets  can  therefore  form  a  complete  set  of  basis  functions  for  any  function  that  has 
finite  energy. 


Fet  f(t)  represent  a  function  that  has  finite  energy.  The  function/ft)  can  be  approximated  by  a  linear 
combination  of  the  wavelets.  It  has  been  shown1  that  all  the  features  of  the  function/ft)  that  are  larger 
than  the  scale  2m°  can  be  approximated  by  a  linear  combination  of  the  translates  of  the  scaling  function 
which  is  derived  from  the  wavelet.  Fet  this  approximation  be  A  fit).  Then  we  have 

fit)  =  A  f{t)  +  dm  {t)  where  d  it)  is  called  the  detail.  The  approximation  A  fit)  can  then  be 

decomposed  into  its  approximation  at  a  higher  scale  and  its  detail.  The  approximations  of  the  function 
at  different  scales  can  be  obtained  as  outlined  above. 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 
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APPROACH 


We  adopted  two  approaches  for  obtaining  the  estimates  of  the  sediment  compressional  wave  speed 
profile  at  different  scales.  In  the  first  approach,  the  pressure  field  as  a  function  of  range  obtained  with  a 
monochromatic  source  is  used  as  data.  Let  p(r,a> )  be  the  measured  pressure  field.  Using  wavelets, 
approximation  to  the  pressure  field  at  different  scales  is  first  obtained.  This  approximated  data  can  then 
be  used  to  infer  the  sediment  acoustic  properties. 

In  the  second  approach,  we  use  the  scaling  function  at  different  scales  to  represent  the  unknown 
compressional  wave  speed  profile.  The  unknown  function  c(z),  which  represents  the  compressional 
wave  speed  in  the  sediments,  can  then  be  approximated  to  scale  m  by  the  relation 

AmC(Z.)  =  Y.an(f>m,n(Z) 

where  (z)  is  the  scaling  function  at  scale  m  and  Amc(z)  is  the  approximation  of  c(z)  to  scale  m. 
Inserting  this  expression  into  the  inversion  algorithm,  we  obtain  the  coefficients  a”'  which  are  then  used 
in  the  equation  above  to  obtain  the  compressional  wave  speed  profile  at  scale  m. 

WORK  COMPLETED 

The  scheme  was  tested  with  synthetic  data.  We  consider  a  range  independent  model  of  the  ocean.  The 
compressional  wave  speed  profile  in  the  sediment  and  the  water  column  are  shown  in  Figure  1.  For  this 
model  and  assuming  a  known  constant  density  in  the  sediments,  the  pressure  field  as  a  function  of  range 
was  obtained  for  a  source  frequency  of  175  Hz.  Approximations  to  the  pressure  field  obtained  using 
Daubechies  wavelet  for  scales  7,  6,  and  5  and  the  corresponding  wavenumber  spectrum  are  shown  in 
Figure  2. 
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Figure  1:  Sediment  compressional  wave  speed  profile  and  sound  speed  profile 
in  water  column.  Depth  of  water  column  =  70  m. 
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Figure  2:  True  field  and  scale  7  and  5  approximations  to  the  field  and  their  spectrum. 


Figure  3:  Sediment  compressional  wave  speed  profile  from  different  approximation  to  the  field. 


We  note  from  Figure  2  that  the  number  of  propagating  modes  for  the  true  field  is  14.  In  the  scale  7 
approximation  to  the  field  there  are  only  four  propagating  modes.  However  the  scale  5  approximation 
of  the  field  is  in  close  agreement  with  the  true  field.  Also  the  wavenumber  spectrum  is  similar  to  that  of 
the  true  field  though  the  modal  amplitudes  of  the  higher  order  modes  are  less.  The  different 
approximations  of  the  field,  therefore,  corresponds  to  spatially  filtering  the  data  with  a  band  pass  filter. 
The  bottom  compressional  wave  speed  profile  corresponding  to  these  scales  are  presented  in  Figure  3. 
We  note  that  at  scale  7,  the  inversion  results  in  estimation  of  the  sediment  compressional  wave  speed 
profile  at  very  shallow  depths.  The  inversion  at  scales  6  and  5  yield  estimates  of  the  compressional 
wave  speed  in  the  deeper  layers.  Similar  analysis  was  performed  using  the  data  acquired  during  the 
Hudson  Canyon  experiment  conducted  in  1988.  The  results  of  this  analysis  are  shown  in  Figure  4. 

In  the  second  approach,  we  obtained  for  the  ocean  model  described  earlier  the  eigenvalues  of  the 
propagating  modes  at  50  Hz,  100  Hz  and  200  Hz  using  KRAKENC2.  The  use  of  modal  eigenvalues  for 
estimating  the  compressional  wave  speed  profile  in  sediments  has  been  well  documented3.  In  our  earlier 
work,  boxcar  functions  were  used  to  represent  the  compressional  wave  speed  profile  in  the  sediments. 

In  this  study  we  use  scaling  function  derived  from  wavelet  to  represent  the  unknown  compressional 
wave  speed  profile.  The  scaling  function  derived  from  Daubechies  wavelet  was  used  in  the  analysis. 

The  results  of  reconstruction  showing  the  results  for  different  scales  are  presented  in  Figure  5. 

RESULTS 

The  approximation  to  the  field  using  wavelets  yields  spatially  filtered  field  in  which  at  the  highest  scale 
only  the  low  order  modes  exist.  This  provides  an  estimate  of  the  acoustic  properties  close  to  the 
sediment  water  interface.  As  move  to  lower  scales,  the  inversion  yield  properties  at  deeper  depths.  In 
the  second  method  where  we  used  multi- frequency  data  the  maximum  resolution  attainable  is  related  to 
the  highest  frequency  of  the  source  and  is  approximately  half  the  wavelength  at  this  frequency. 

IMPACT/APPLICATIONS 

Multi-resolution  method  permits  us  to  obtain  the  best  possible  estimate  of  the  bottom  properties  that  the 
data  will  support.  It  can  be  used  to  look  selectively  at  different  regions  and  obtain  estimates  with  high 
resolution  close  to  the  water  sediment  interface  and  estimates  with  coarser  resolution  at  deeper  depths. 

RELATED  PROJECTS 

During  the  period  of  the  contract,  additional  work  on  inversion  for  density  profile  from  modal 
eigenvalues  was  completed.  Also  a  comparison  of  local  and  global  algorithms  for  inversion  of  bottom 
acoustic  properties  was  performed. 
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Figure  4:  Measured  field  (175  Hz)  and  field 
predicted  by  model  from  scale  5 
approximation  (Top).  Model  from  scale  5 
approximation  of  field  and  model  obtained 
from  multi-frequency  inversion 4  (Bottom). 


Figure  5:  Sediment  compressional  wave 
speed  profile  at  different  scales  from  multi¬ 
frequency  data. 
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